NASA Technical Memorandum 107206 


Predicting the Stability of a Compressible 
Periodic Parallel Jet Flow 


Jeffrey H. Miles 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

1996 Fluids Engineering Conference 

sponsored by the American Society of Mechanical Engineers 
San Diego, California, July 7-11, 1996 



National Aeronautics and 
Space Administration 



Predicting the 
Stability of a Compressible 
Periodic Parallel Jet Flow 

Jeffrey Hilton Miles 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, OH 44135 


Abstract 

It is known that mixing enhancement in compressible 
free shear layer flows with high convective Mach num- 
bers is difficult. One design strategy to get around this 
is to use multiple nozzles. Extrapolating this design 
concept in a one dimenional manner, one arrives at an 
array of parallel rectangular nozzles where the smaller 
dimension is w and the longer dimension 6 is taken to 
be infinite. In this paper, the feasibility of predicting 
the stability of this type of compressible periodic par- 
allel jet flow is discussed. The problem is treated using 
Floquet-Bloch theory. Numerical solutions to this eigen- 
value problem are presented. For the case presented, 
the interjet spacing, s, was selected so that s/w = 2.23. 
Typical plots of the eigenvalue and stability curves are 
presented. Results obtained for a range of convective 
Mach numbers from 3 to 5 show growth rates u>i = *2* 
range from 0.25 to 0.29. These results indicate that 
coherent two-dimensional structures can occur without 
difficulty in multiple parallel periodic jet nozzles and 
that shear layer mixing should occur with this type of 
nozzle design. 

1 Introduction 

Interest in proving the economic and environmental fea- 
sibility of a high-speed civil transport has stimulated 
again studies in supersonic nozzle design. The primary 
process which must be understood and controlled to 
meet design objectives involves supersonic mixing. Im- 
proving the efficiency of supersonic combustors is an- 
other area where mixing enhancement at supersonic 
speeds must be understood and controlled to meet de- 
sign objectives. 

A reduction in mixing and a growth rates of su- 
personic shear layers with increasing Mach number has 
been demonstrated experimentally by many investiga- 


tors (Refs. [1], [2], [3] ). In addition, linear stability 
analysis shows results that are similar to the experi- 
mental studies (Refs. [4], [5], [6] ). However, these 
linear stability investigations have been of unbounded 
flows that are uniform at infinity. No attention has 
been given to supersonic flows periodic in the direction 
normal to the stream. In this paper, the feasibility of 
predicting the stability of a compressible periodic par- 
allel jet flow is discussed. While this problem has not 
been dealt with analytically, one design strategy to get 
around the decrease in mixing with increasing convec- 
tive Mach number is to use multiple nozzles. Recent 
flow studies using multiple rectangular nozzles are pre- 
sented in Refs. [7] and [8]. Extrapolating this design 
concept in a one dimensional manner, one arrives at an 
array of parallel rectangular nozzles where the smaller 
dimension is w and the longer dimension b is taken to 
be infinite. In addition, one starting point for a compu- 
tational fluid dynamic study of multiple nozzles might 
be a study of a compressible periodic parallel jet flow. 
Consequently, the work discussed might be of interest 
to both the experimental and computational nozzle de- 
sign community. Furthermore, this type of study might 
also be useful in understanding screech from multiple 
nozzles since it provides information on the flow distur- 
bances leaving the nozzle. 

In this paper the linear stability of a spatially peri- 
odic supersonic jet flow having velocity profile charac- 
terized by an abrupt rise and a sharp fall is discussed. 
In contrast to the uniform velocity distribution of free 
flows such as jets, wakes, and free shear layers in the far 
field, the periodic flow repeats exactly its velocity pro- 
file endlessly. The periodic flow is, therefore, expected 
to have different stability characteristics than those of 
free flows. In particular, the variation growth rates with 
Mach number needs to be examined. In this paper, the 
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inter jet spacing, 5, was selected so that s/w — 2.23 and 
the ratio of the periodic nozzle wavelength, X — s -\-w, 
to w was 3.23. Note that w is the smaller dimension of 
the rectangular nozzle. Results obtained for a range of 
convective Mach numbers from 3 to 5 will be presented. 


where the primes denote differentiation with respect 

to y 
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2 Formulation of the problem 

Let (U(y), 0, 0) be the velocity of a steady plane-parallel 
flow, where the x-axis is in the direction of the flow and 

U(y)=U + —^—h(y) 

where U\ is the velocity outside the jet, U 2 is the 
mean centerline jet velocity, U — , A U — U^ — U 1 , 

and h(y) is the velocity profile function which varies 
from -1 to 1. 

The flow field is perturbed by introducing wave dis- 
turbances in the velocity and pressure with amplitudes 
that are a function of y. Thus, 

(i *,v,w,p) 

= («(y),v(y),w(y),p($))exp |i (fci + /f-wr)j . 

Where 


In this paper, the velocity profile function, h(y), is 
periodic such that 

A(y +2ir) = h(y). 

The velocity profile h(y) is not any exact solution 
of the Navier-Stokes equation, but it can be considered 
as a simple model of some real periodic flow. 

The velocity profile h(y) discussed herein is given 

by 

My)=l-2/(y) 
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and from Crocco’s Equation [9] 
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By definition k is real positive number that rep- 
resents the wavenumber in the x-direction, £ is the 
wavenumber in the z-direction, Cr is the relative phase 
velocity, and u>i = ^ is the amplification rate of the 
disturbance. 

From the equations of motion if nonlinear and 
viscous terms are neglected one can obtain an equation 
for the y-component of the perturbation velocity as 
follows: 


V 


T* A ' 

»'<T + T> 

T A 


( 2 . 1 ) 


h" 


( h-c ) 


+ Aik 


T' A' h' 

( f + A ’(h-c) 


= 0 


where the function /(y) is given by 


/(y) = 


1 + ( sinh (^rr)) 


18 


TI = A(-l + 7), 

A = 1.5 and y goes from 0.0 to 2n. The profile 
function f(y) is adapted from an equation used by 
Monkewitz in Ref. ( [10] ) in a study of the absolute and 
convective instability of two-dimensional wakes. Only 
two-dimensional disturbances will be considered. A 
schematic of the nozzle geometry is shown in Figure 
1. The velocity profile is shown in Figure 2. 
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3 Floquet-Bloch theory 

Since the basic flow velocity profile , /(y), is periodic, 
equation (2.1) is an example of a Floquet-Bloch prob- 
lem. The mathematics of solving Floquet-Bloch type 
problems is discussed in Refs. [11], [12], and [13]. Ap- 
plications to solid state physics are discussed in Refs. 
[14], [15], and [16]. Applications to spatially periodic 
flow is discussed in Refs. [17], [18], [19], [20], and [21]. 

The paper by Beaumont (1981) Ref. [19] and the 
description of the Floquet-Bloch theorem by Hochstadt 
(1963) in Ref. [12] were particularly useful in guiding this 
research. 

A survey of the spatially periodic flow literature is 
presented by K. Gotoh and M.Y. Yamada in Ref. [22] . 

The second order differential equation can be de- 
scribed by a system of first order differential equations. 
Let 

v = X\ 

V* = X2 

so that Eq. 2.1 can be rewritten as the system 

(3- 2 ) X ' = ( D C ) X 

where 


and 


|*(2jt) - pl| 

_ 4> i(2tt) - ^ <^2(27r) 

fa' (2 k) ~ A* 

= /* 2 — (<£i(2*r) + <^2 , (27r))^ 

+(^ 1 (27r)0 2 '(2 5 r) - 

^ 2 ( 2 *)^!/ (2?r)) 

(3.3) = n 2 - (4>i(2*) + <f>2( 2*))m + 1 

Since 

<h(2*)^(2*) - <M2*)^i'(2tt) = |*(2ir)| = |«(0)| = 1 

The independent solutions of equation (3.2) have 
the form 

<t> = X(y) ex P( ~|~ y) = x (v) exp( r v) 

The parameter T specifies the period of the eigen- 
function <f>. If T is real the eigenfunction grows or de- 
cays at infinity. Consequently, only imaginary values 
of T are acceptable. Thus the eigenfunction oscillates 
in space and is called a continuous mode. The distur- 
bance with Ti = 1/n, where n is a nonzero integer, has 
a period 2n7r. One with T* = 0 has the same period 2ir 
as the main flow, while an irrational value of r\ means 
the disturbance is aperiodic. Note that the parameter 
T does not appear in the flow equation, but is due to 
the Floquet-Bloch theorem. 

Solutions of 3.2 are thus of the form 


D = 


ft" 

{h-c) 


- X 

+ Aik - (^r + 
T 



h' 

(h-c) 


If*(y) is a fundamental matrix solution of equation 
( (3.2 ) ) such that 


*(0) = I 

where I is the identity matrix, then from the 
Floquet-Bloch theorem 


*i(v + 2tt) = Pi*i(y) 

Xi(y+2n) = M 2 -X 2 (if) 

where m and y 2 represent the zeros of (3.3), pro- 
vided they are distinct. 

In general, these solutions will not be periodic. 
Conditions for periodic solutions can be found as 
follows 

Let fxi = e i$i and /12 = e~ t0i . 

Then from equation (3.3) 


$(y +2 tt) = 4>(y)$(27r) 


We now introduce two solutions of equation ( 3.2 ) 
with initial values at y = 0.0 
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Next we seek the eigenvalues of <fr(27r) 


cos {0 % ) = <j> i(27r) + <j> 2 / (27r) — 6/2 

Consequently, for a solution to be periodic 6 must 
be real and \6\ smaller than 2. 

The constants [x are termed the characteristic mul- 
tipliers of the Floquet-Bloch system (3.2) and the cor- 
responding characteristic exponents are determined by 
the relation T = T r + = f£ + * |** 
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Table I Eigenvalues at maximum growth rate 


(d = 0.2 A = 1.5 s/h = 2.23) 
I i u I I r. 


1712 

k 


r< 

Cr 

3 

2.855 

0.2855 

2.158e-01 

-0.2195 

3.2 

2.735 

0.2735 

1.359e-01 

-0.2005 

3.4 

2.495 

0.2495 

9.450e-02 

0.4963 

4. 

2.470 

0.2470 

7.632e-02 

0.5262 

5 

2.675 

0.2675 

4010e-03 

0.6367 


Table II Cray YMP CPU 

computation time, hrs. (A = 1.5 s/h == 2.23) 


1712 

<* = 0.1 

P 

II 

O 

to 

to 

© 

II 

2 

Total 

3 

20.863 

12.768 

3.049 

36.68 

3.2 

20.654 

11.697 

6.913 

39.265 

3.4 

19.236 

12.724 

5.571 

37.531 

4 

26.147 

11.847 

6.525 

44.519 

5 

25.555 

13.17 

5.212 

43.937 

Total 




201.932 


4 Numerical calculations 

In obtaining the results presented here equation 3.2 
was integrated from y = 0 to y = 2 tt using a standard 
Runge-Kutta procedure (IMSL math library routine 
IVPRK) for fixed values of c* = 0.1, 0.2, and 0.3 for 
m2 values of 3, 3.2, 3.4, 4 and 5. 

To investigate stability, iteration was used to vary 
the real value of the phase velocity, Cr. At each iteration, 
the eigenvalues, p, of the matrix $(27 t) were then 
found and the corresponding characteristic multipliers 
r calculated. The iteration was successful if Si = 0 so 
that T r = 0. Also at each iteration a check was made 
that det $(2n) = 1 . 

Typical eigenvalue plots of the phase velocity, Cr 
verses growth rate u>i and stability plots of the charac- 
teristic multipliers Tt verses are presented in Figures 
3-8 . 

5 Results 

Results for m2 = 3 are given for c, = 0.1, 0.2, and 0.3 in 
Figures 3-5. The (a) part of the figure shows eigenvalue 
plots of the phase velocity, dr verses growth rate u ;» while 
the (b) part of the figure shows stability plots of the 
characteristic multipliers I\ verses ui. Results using 
the same format are presented for m2 = 4.0 in Figures 
6 - 8 . 

Table I shows the values at the growth rate maxi- 
mums for the cases investigated. Results obtained for 
a range of convective Mach numbers from 3 to 5 show 
growth rates cj* = range from 0.2 to 0.3. 

Table II. shows the CRAY YMP CPU compu- 
tational time needed to search the phase speed, c, 


wavenumber k space in order to see if solutions exist. 
The wave number spacing used was 0.005. The total 
computation time was slightly more than 200 hours. 
The results were slowly accumulated over a period of 
seven months. 

6 Conclusions 

A numerical method for performing the stability analy- 
sis of a compressible periodic parallel jet flow has been 
developed and tested. Obtaining a set of solutions is a 
tedious expensive process. However, the resulting solu- 
tions provide detailed information on the growth, phase 
speed and periodicity of small disturbances propagating 
with the flow near the nozzle exit. 

These results indicate that coherent structures can 
occur without difficulty in multiple parallel periodic 
jet nozzles and that shear layer mixing should occur 
with this type of nozzle design. Using these results, 
it should be possible to determine which nozzle design 
yields solutions with high growth rates and to pick out 
likely screech frequencies due to feedback from shock 
cells to the nozzle. 
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Figure 8. — (a) Eigenvalue c r verses growth rate ojj = kc/2 (Cj = 0.3, m 2 = 4, s/w = 2.23). (b) r s 
growth rate toj = kc/2 (Cj = 0.3, m 2 = 4, s/w = 2.23). 
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